function equation = get_equation_of_reaction(RKEGGID)
    reactions = evalin('base', 'reactions'); 
    for i = 1:length(reactions)
       if strcmp(reactions(i).ENTRY, RKEGGID)
           equation = reactions(i).EQUATION;
           break;
       end
    end
end